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Abstract 

Through introducing a notion of renormahzation of particle-number 
density, a simple perturbation scheme of nonequilibrium quantum-field 
^ ■ theory is proposed. In terms of the renormalized particle-distribution 

d \ functions, which characterize the system, the structure of the scheme 

(and then also the structure of amplitudes and reaction rates) are the 
same as in the equilibrium thermal field theory. Then, as an obvious 
consequence, the amplitudes and reaction rates computed in this scheme 
are free from pinch singularities due to multiple products of (5-functions, 
which inevitably present in traditional perturbation scheme. 
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Since late fifties, efforts liave been made to incorporate quantum field tfieory witli 
nonequilibrium statistical mechanics |Q 0, ^. In connection with the description of 
quark-gluon plasma, which is expected to be produced in heavy-ion collisions and to 
have existed in the early universe, much work has recently been devoted to this issue 
(see, e.g., § g § 0, g |). 

The standard framework of nonequilibrium theory is formulated by employing the 
closed-time path in a complex-time plane. A few years ago, Altherr and Deibert 
have pointed out that, when calculating amplitudes or physical quantities in pertur- 
bative closed-time-path formalism, the "pathological terms" p, Q necessarily appear. 
More precisely, such terms appear in self-energy-parts-inserted two-point functions 
in association with pinch singularities coming from multiple products of 5-functions. 
These terms do not cancel each other unless the particle distributions are those for 
a system in thermal and chemical equilibrium. Since then it is shown that classes of 
such pathological terms can be resummed (see also [0). An application of this re- 
sult to the rate of hard-photon production from a nonequilibrium quark-gluon plasma 
system is made in [Q. 

In this letter, introducing a notion of renormalization of number density, we pro- 
pose a simple perturbation scheme of quantum field theory for computing the rates 
of reactions taking place in a quasiuniform system near equilibrium or nonequilib- 
rium quasistationary system. The scheme is the same in structure as the equilibrium 
thermal field theory (ETFT), provided that the equilibrium number density n{pQ) 
in ETFT is replaced by the renormalized nonequilibrium number density N{P) that 
characterizes the system under consideration. Then, in contrast to the canonical per- 
turbation scheme, no "pathological terms" mentioned above appear in amplitudes 
and reaction rates evaluated within the scheme. 

We employ the closed-time-path formalism p|, |^, In this formalism, propaga- 
tors, vertices and self-energy parts enjoy 2x2 matrix structure. We use 13 to denote 
the 2x2 matrix whose (ij) component is Bij. A vertex matrix V has simple struc- 
ture, Vi2 = V21 = and Vu = —V22, with Vu the vertex factor in vacuum theory. 
Let A(a;i,X2) and S(xi,X2) be, in respective order, the bare propagator and the self- 
energy part in a configuration space. Following p[, we assume that reactions taking 
place in a system under consideration are "described" by the relative or microscopic 
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coordinates x = xi — X2- Then, we make Fourier transforms with respect to x'^, 
A{P,X) and S(P, X), where X = [xi +X2)/2 is the center-of-mass or macroscopic 
coordinates. It is assumed that the dependence of A{P,X) and S(P, X) upon X^^ 
is weak. More precisely, over the macroscopic space-time region of the system, where 
a microscopic or elementary reaction takes place, X^-dependence of A(P, X) and 
E(P, X) may be ignored. Then, in calculating the reaction rate, we can use A(P, X) 
and S(P, X) with fixed X^. This is the first stage of the theoretical analysis. Mi- 
croscopic reactions cause changes in the number densities of (quasi)particles, through 
which the density matrix changes with macroscopic space-time X^. Dealing with 
this is the subject of the second stage, where (weak) X^'-dependence of A(P, X) and 
S(P, X) are explicitly taken into account. In this letter, as in [0, |§, 0, we concentrate 
our concern on the first stage and drop the argument X throughout. 

The matrix elements of A(P) and S(P) enjoy various properties 0]. Among 
those we shall use ReE^ = PeSai = 0, S22 = -^u, E?,j=i(-)'+^'Aii(^) = and 
i:S,=iS.,(P) = 0. 

For simplicity of presentation, we take a massless, self-interacting, complex-scalar 
field theory with a conserved charge. A(P) may be written as 

A(P) = M(P)A^(P)M(P), (1) 

whereQ 

A^(P) = diag(A^(P), -A^(P)) (A^(P) = l/(p2 + ^e)) , (2) 

e{-po)+n{P) 

'-M+m Jl + n{P) 



1 + n(P) 

M{P) = I ' .„..„.p. I . (3) 



Here n{P) = n{pQ, p) with p = |p| is the number-density function that characterizes 
the ensemble of the systems. In the case where the "local" temperature T = (3~^ 
and the chemical potential /x, being conjugate to the charge, are defined, n{P) = 
]^/|^g/3(|po|-e(po)M) p and ^ the functions of macroscopic space-time coordinates. 

Here we recall that, in ETFT, the self-energy-part matrix S(P) takes |^ the form 

S(P) = M-\P) tpiP) M-\P) tp{P) = diag (S^(P) - S^(P)) . (4) 



t In calculating an amplitude in ETFT, a common practice |Q is to keep e finite throughout and 
at the end of calculation the limit e ^ 0+ is taken. We follow this procedure. 
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One can easily see from Eqs. ([^), (0) and (^) that an / (> 1) self-energy-parts- 
inserted two-point function Gi{P) = A(P)[S(P) A(P)]' does not include functions 



(p2 ^ i^yk^p2 _ 



-{l+l-k) 



[1 < k < I) but includes functions (P^ ± ie) In 



the limit e 0+, the former functions develop pinch singularity while the latter 
functions turn out to the well-defined distributions, and thus G/(P) is free from pinch 
singularity. In the present nonequilibrium case, Eq. (^) does not hold and Gi{P) 



It 



It 



I '^^^^ . This is essentially 



includes "pinch-singular" functions (P^ 
what has been observed in 0. 

Let us construct the renormalization theory by introducing a renormalized or 
"physical" number density N{P), 



N{P) = n{P) + 5n{P) , 
with which the bare (renormalized) propagator reads (cf. Eqs. (|Tp - (^) 

A(")(P) = A(P) 
m(")(p) = M(P) 



(5) 



n(P) N{P) 



(6) 



n(P) -* N{P) 

In order to compensate the difference between the renormalized A^*") and the original 
"bare" A, 6A{P) = A^^^{P) — A(P), we should introduce the counter term in the 
interaction Lagrangian, 

1 /• ^ 2 

Lc=-- / n (^"^^0 ^ 0*(a;i) [A"^(xi, Xs) 5A(X2,X3) A"^(X3, X4)j ..0j(x4) , (7) 

^ e=i i,j=i 
where (pj and 0* (j = 1, 2) stand for the type-j fields and 



A-\x,y) = 
SA{x,y) = 



d^P 
d^P 
d^P 



g-iP.(x-s/)^-l^p^ 
^-iP-{x-y) ^^^p^ 



-iP-{x-y) 



Here Aj 



(27r)4 

1 (i, j = 1, 2) for all i and j and^ 5e(-P^ 



27Ti 6n{P) 6,{P'^)A 



./[7r{(P 



2\2 



The perturbative calculation goes with propag ator A('^)(P) in Eq. (D and ' ver- 
tices" coming from Li^t + L^.., where Li^t is the original interaction Lagrangian. Con- 
forming to the ultra-violet renormalization theory, one can say that n{P) is the bare 
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number-density function and N{P) is the renormalized or physical number-density 
function. The latter describes the physical system under consideration. The self- 
energy part S*^'')(P) obtained from the above perturbation scheme takes the form^ 

SM(P) = t{P) - A-\P) 6A{P) A-\P) . (9) 

We are now in a position to set up the renormalization condition: The self- energy 
part S('')(P) zs the same in structure as Eq. ^ in ETFT, i.e., 

^^(P) = [mW(P)]-i ±P{P) [M('-)(P)]-i (10) 

where 

S?(P) = diag (e?(P), -(S?(P)) j . (11) 

The second term on the right-hand side of Eq. (P) does not contribute to {TiP)ii 
and (S?)22 and the condition [{S?(P)}22]* = -{S?(P)}n is automatically met; 

sP(p) = Sn(P) + e{po) Si2(P) + e{-po) S2i(P) . (12) 

From Eq. (||) with Eq. (^ and Eqs. (H) and (^), we see that the condition 
{eS:)(P)}i2 = {S?(P)}2i = yields 

6n{P)6,{P') = -^p^^^52^[^^(Po){(l + iV(P))Si2(P)-iV(P)S2i(P)} 
+9{~po) {(1 + iV(P))S2i(P) - iV(P)) Si2(P)}] . (13) 

Now let us introduce 

rj-*)(p) = eipo)[{i + NiP)}r,{P)-N{P)UP)] 

+ei-po)[{l + NiP)}T,{P) - N{P)T,{P)] , (14) 

where 

r,(P) = r,(P) = -^Si2(P), 

r,(P) = r,(P) = -^S2i(P). (15) 

^ Suppose that we are analyzing in the nth order of perturbation series. The second term on 
the right-hand side of Eq. (|^) comes from with the nth order counter term, being proportional 
to (5n(")(P). The nth order E(P) in Eq. (|) involves Sn^^^P) with 1 < j < n - 1. 
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On the mass shell po = p (po = —p), r^"'^*-'(P) is the net production rate of a particle 
(antiparticle) with momentum p (— p) [0]. Then T^^^\P) in Eq. is an off-shell 
extension of the net production rate. Substitution of Eq. with Eq. (|T5|) into Eq. 
(0) yields 

M^)^e(n = ^(p^y^r(-)(P). (16) 

It is to be noted in passing that, in ETFT, T^^^\P) vanishes (see, e.g., p[). 
In the limit e 0"*", 6n{P) in Eq. ([T6|) is singular due to the pinch singularity at 
|po| = P- Equation ( [T^ ) may be "solved" for 6n{P) in the form, 

/+00 //t-r, 
A^ko, P) A^iko, p) (17) 
-oo /vr 

= ^r("^*)(P), (18) 

where e ~ e/2p. 

It is obvious from the above construction that the renormalization theory is ob- 
tained from the ETFT simply by substituting N{P) for the equilibrium distribution 
function neipo) (= l/[e^{|pol-^(po)^) - 1]). Then, the fact in ETFT that a self-energy- 
parts- inserted propagator is free from pinch singularities (cf. above after Eq. (^) 
is transmitted as it is to the case of renormalization theory. In fact, from Eq. (^) 
with Eqs. (g) and (gj, and Eqs. (|TUp and (|n|), we see that M'-\P)[t^'-\P) M'-\P)f 
(/ > 1) includes only well-defined function (P^ ± ze)"^'^^-*. 

We now demonstrate that, as far as the singular pieces are concerned, n(P) on 
the mass shell po = p [po = —p] coincides with the Heisenberg particle [antiparticle] 
number density n^{p) [n^ {p)]. n^{p) reads, with obvious notation, 

(at(p, t)a(p', t)) = 5(p-p')nf(p) 

1 f d^xd^y ■fp,^_p,. d d 



= yo = t 

(19) 



where P = (p, p) and P' = (p', p') . For the case of antiparticle number density, we 
have Eq. (||) with 0* 4-^ 0. Note that (0*(x)0(?/)) [(0(x) 0*(?/))] is the (12) [(21)] 
component 0, H, Q of the full propagator iGi2{y — x) [iG2i{x — y)], we obtain 

na/biP) = YpJ^^ 27 ^ ^G'i2/2i(A;o, p) • (20) 
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We compute n^/hip) up to the contribution from one self-energy-part-inserted 
6*12/21(^0) p)- We first employ the perturbation scheme in the renormalization theory 
defined above. As mentioned above, computation goes just as in the case of ETFT. 
The resultant n^/f,{p) reads 

na/bip) = N{±p, p) + 5N{±p, p) . (21) 

Here comes from the lowest-order contribution to G12/21 and 6N comes from the 
one self-energy-part-inserted Gi2/2i- As in the case of ETFT, 6N is a well-defined at 
most finite functional of A^. 

Now we turn to compute n^/^ip) in traditional scheme, without introducing the 
renormalization counter term. The lowest-order contribution to 6*12/21 yields n^/h{p) 
= n{±p, p), where n is the "starting" number density function in Eq. (^. Then, we 
write 

n^/,{p)=n{±p,p) + Sn^/,{p). (22) 
The contribution from one self-energy-part-inserted 6*12/21 to 5n^i^{p) reads, 

^(fco±p)^[A^Apf("^*)(A;o,p) 
-{6^ + n) {{/^pftp - (A^)2S^}] , (23) 

where T^p = Sn + 6'+Si2 + ^-S2i, d± = 6{±ko), n = n{ko, p), Ap = Apiko, p) and 
= ^F{ko, p). T^^^^\ko,p) in Eq. (|23| ) is the traditional-scheme counterpart of 
rJ,"^*)(P) in Eq. ([1|) with Eq. (|15D. In the limit e ^ 0+, the term with A^Ai. 
in Eq. (|23| ) diverges due to pinch singularity at |fco| = P- On the other hand, the 
contributions from (A^?)^ and {A*py, being well-defined functions, are at most finite. 
We deal only with the singular contribution coming from |/co| = p, 



(p) ^ 2p' *)(±p, p) ^A^(A;o, p) Apiko, p) , (24) 

J-00 ZTT 



where '~' is used to denote an approximation that is valid for keeping the singular 
contribution. In obtaining Eq. (0), use has been made of the fact that A*p Ap is even 
function of k^. 5n{P) on the mass shell Pq = ±p, Eq. (|T^), and the singular piece of 
5n^l^{p) in Eq. (|2^ ) is the same in form. It should be noted, however, that rj,""^*-* 7^ 
-p(net)_ rewrite r^^*^*) in terms of A^ by using n(±p, p) = N(±p, p) + {ri(±p, p) — 
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N{±p, p)}. The philosophy of renormahzation tells us that the difference n — N only 
affects the higher order contribution to n^/b{p)- Thus, to the accuracy of one self- 
energy-part insertion, we have 5n^i^{p) ~ 5n(±p, p). Then, from Eqs. (|2T|) and (|2^) , 
we find 



^f/bb) ^N[±p,p) ~ n{±p,p) + Sn^/t,{p) 



n{±p, p) + 6n{±p, p) , (25) 



where use has been made of the fact that 6N , Eq. (21), is not singular. The relation 
(P5|) is in accord with Eq. (^ on the mass shell po = ±p. 

Thus, we have learned that, as far as the singular pieces are concerned, the renor- 
malized number density N{P) (Eq. (^) on the mass shell po = ±p coincides with 
the Heisenberg number density n^/h{p)- Then we may regard N{P) as an off-shell 
extension of the Heisenberg number density. 

In passing it is worth making a comment. The first computation leading to Eq. 



(pT|) is based on the theory, in which the free Lagrangian is written in terms of physical 
or renormalized number density N{P). On the other hand, the second computation is 
based on the theory, in which the free Lagrangian is written in terms of original bare 
number density n{P). As in the ultra-violet renormahzation theory in quantum field 
theories, both approaches are equivalent. In the first approach, the renormahzation is 
"done" at the beginning by introducing the counter Lagrangian Lc. While the second 
approach starts with the bare Lagrangian and the renormahzation is "done" at the 
end. 

Now let us turn to analyze the physical meaning of the singular part of 5n^^(p), 
Eq. (p^ . Introducing the Fourier transform of A.p{kQ, p), Eq. (^, 

-°° dh 



Ap{t,p) ^ [ ^^e-^''o'Ap{ko,p) 

J -oo ZTC 

d{t)e-'^^-''^' + d{-t)e'^^-''^'] , (26) 



-oo 

i 



2p 



we obtain 



/+00 Wi^ r+QO 
AUko, P) Ap{ko, p)= dt AU-t, p) Ap{t, p) . (27) 
-oo Z7[ J —oo 

Here let us recall that, in standard Lippmann-Schwinger formalism, a singular func- 
tion 27c[6{Ef — Ei)]"^ appears in a transition probability of some reaction, where Ei 
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(Ef) is the energy of the initial (final) state of the reaction. We follow the standard 
argument for interpreting this function as (tj — ti)6{Ef — Ei), where tf — ti is the time 
interval during which the interaction acts. Applying this to the present case, we have 



dtAU-t,p)AF{t,p) = 



2i 2i 



{2i 



Substituting this into Eq. (pTj) and then into Eq. (El), we obtain 



5n%,{p)^^-^^^r\±P.P)- (29) 

It is to be noted in passing that if we take the limit ti — oo and tf ^ +00 in Eq. 
(p8|), we have l/4p^e ~ l/2pe. Substitution of this into Eq. (p^ reproduces Eq. ( |18|) 
on the mass shell po = ±p. 

It is interesting to note that Eq. (p9|) is half of the net production probability 
during the time interval tf — ti of the reaction. What the result (pOf ) tells us is the 
following. Since F^"*^*) ^ 0, the number density Ua/bip, t) changes with microscopic 
time t. At the initial time ti (~ —00) of some reaction na/b{p, ti) = n{±p, p), Eq. 
(I), and at the final time tf (~ +00) ha/b{p,tf) = ha/b{p, ti) + (t/ - ti)r(,"^*)(±p, p) ~ 
n{±p, p) + 26na/b{p). Then, 5nf/^(p) in Eq. (||) is 6n^/^{p) = na/b{p, {tf + ti)/2) - 
na/b{p,ti), i.e., N{±p,p) ~ n{±p,p) + 5n^n^{p), Eq. (|25|), is the number density 
'^a/bip^t) at the middle or average time of the initial and final times of the reaction, 
t={tf + ti)/2. 

Here we briefly mention how to treat the gauge field and massless fermion field. 
In the case of gauge field like gluon in QCD, we employ the Coulomb gauge or the 
Landshoff-Rebhan's variant [jTU[ of the covariant gauge. The r;,(P)-dependent part 
of A(P) and S(P) are decomposed into two sectors, the electric sector (e) and the 
magnetic sector (m). For each sector r (= e, m), we introduce the Bogoliubov matrix, 
which is given by Eq. (^ with n{P) n^'^\P). For massless fermion field, we 
decompose A(P) and t{P) into two sectors, A(P) = X;r=± ^rA(^)(P) and S(P) = 
E.=±P.S^"n^)> where = (7o-rp-7)/2. In particular, A^(P) = E.=±^rA?(P) 
with AP(P) = diag[l/{po(l + «e)-rp},-l/{po(l-«e)-rp}] (cf. Eq. (D). For each 
sector r (= ±), we introduce the Bogoliubov matrix, Eq. (^, with n{P) — > —rS'^\P). 
With this preliminaries, the analysis goes as in the case of complex-scalar field. 
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Finally we mention two related work. The first one is the so-called nonequilibrium 
thermo field dynamics This is a "single-time" formalism, without distinction be- 
tween the microscopic- and macroscopic-times. The time representation, rather than 
the po representation, is employed and through carrying out renormalization of the 
number density a nonrelativistic Schrodinger field theory is analyzed. Since the time 
representation is used, the structure of the (Schrodinger) propagators in a pq space 
may not be seen directly. The second one starts with performing the approximate 
resummation of the absorptive parts of the self-energy parts in relativistic scalar field 
theories. Then, the perturbation scheme is so constructed that the renormalized free 
Lagrangian Lq"^ yields the approximately resummed propagator. The counter La- 
grangian (the difference between the original free Lagrangian and Lq'^) is local. By 
contrast, the counter Lagrangian (^ in our renormalization scheme is nonlocal as in 
the hard-thermal-loop resummation scheme 0, 11] in ETFT. 
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